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Abstract. We establish existence and pointwise estimates of fundamental so- 
lutions and Green's matrices for divergence form, second order strongly elliptic 
systems in a domain (i C K. n ,n > 3, under the assumption that solutions of 
the system satisfy De Giorgi-Nash type local Holder continuity estimates. In 
particular, our results apply to perturbations of diagonal systems, and thus 
especially to complex perturbations of a single real equation. 



1. Introduction 

In this article, we study Green's functions (or Green's matrices) of second order, 
strongly elliptic systems of divergence type in a domain f2 C M" with n > 3. 
In particular, we treat the Green matrix in the entire space, usually called the 
fundamental solution. We shall prove that if a given elliptic system has the property 
that all weak solutions of the system are locally Holder continuous, then it has the 
Green's matrix in fi. (For example, if coefficients of the system belong to the 
space of VMO introduced by Sarason [19] . then it will enjoy such a property). For 
such elliptic systems, we study standard properties of the Green's matrix including 
pointwise bounds, LP and weak LP estimates for Green's matrix and its derivatives, 
etc. 

For the scalar case, i.e., a single elliptic equation, the existence and properties 
of Green's function was studied by Littman, Stampacchia, and Weinberger |14j 
and Griiter and Widman [TT]. In this article, we follow the approach of Griiter 
and Widman in constructing Green's matrix. The main technical difficulties arise 
from lack of Harnack type inequalities and the maximum principle for the systems. 
The key observation on which this article is based is that even in the scalar case, 
one can get around Moser's Harnack inequality 18J or maximum principle but 
instead rely solely on De Giorgi-Nash type oscillation estimates [5] in constructing 
and studying properties of Green's functions. From this point of view, this article 
provides a unified approach in studying Green's function for both scalar and systems 
of equations. We should point out that there has been some study of Green's matrix 
for systems with continuous coefficients, notably by Fuchs [7] and Dolzmann-Muller 
[5J. Our existence results and interior estimates of Green's function will include 
theirs, since as is well known, weak solutions of systems with uniformly continuous 
(or VMO) coefficients enjoy local Holder estimates. On the other hand, we have 
not attempted to replicate their boundary estimates, which depend in particular on 
having a C 1 boundary. Our method does not require boundedness of the domain 
nor regularity of the boundary in constructing Green's matrices, while the methods 
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of Fuchs [7] and Dolzmann-Muller [6] require both boundedness and regularity of 
the domain at the very beginning. We note that a scalar elliptic equation with 
complex coefficients can be identified as an elliptic system with real coefficients 
satisfying a special structure, and thus our results apply in particular to complex 
perturbations of a scalar real equation. In the complex coefficients setting, the main 
results of Section [3] in our paper can be also obtained by following the method of 
Auscher [2] . The estimates of the present paper will be applied to the development 
of the layer potential method for equations with complex coefficients in [T]. 

The organization of this paper is as follows. In Section^ we define the prop- 
erty (H), which is essentially equivalent to De Giorgi's oscillation estimates in the 
scalar case, and introduce a function space Y ' (fi) which substitutes W ' (Cl) in 
constructing Green's functions; they are identical if CI is bounded but in general, 
Yq' 2 (CI) is a larger space and is more suitable for our purpose. In Section [3l we 
study Green's functions defined in the entire space, which are usually referred to 
as the fundamental solutions. The main result is that for a system whose coeffi- 
cients are close to those of a diagonal system, the fundamental solution behaves 
very much like that of a single equation. In Section^ we study Green's matrices in 
general domains, including unbounded ones. We also study the boundary behav- 
ior of Green's matrices when the boundary of domain satisfies a measure theoretic 
exterior cone condition, called the condition (S). We prove in particular that if the 
coefficients of the system are close to those of a diagonal system, then again the 
boundary behavior of its Green's function is much like that of a single equation. 
In section [5j we discuss the Green's matrices of the strongly elliptic systems with 
VMO coefficients. By following the same techniques already developed in the pre- 
vious two sections, we construct the Green's matrix in general domains including 
the entire space. One subtle difference is that in this VMO coefficients case, one 
should play with a localized version of property (H) since basically, the regularity of 
weak solutions of the systems with VMO coefficients is inherited from the systems 
with constant coefficients when the scale is made small enough. Therefore, all the 
estimates for the Green's matrix stated in this section are only meaningful near a 
pole. 

Finally, we would like to mention that when n = 2, the method used in this 
article breaks down in several places and for that reason we plan to treat the two 
dimensional case in a separate paper. 

2. Preliminaries 

2.1. Strongly elliptic systems. Throughout this article, the summation con- 
vention over repeated indices shall be assumed. Let I be a second order elliptic 
operator of divergence type acting on vector valued functions u = (it 1 , . . . , u N ) T 
defined on M™ (n > 3) in the following way: 

(2.1) Lu = -D a (A afj D fj u), 

where A a/3 = A a ^(x) (a, j3 = 1, . . . , n) are N by N matrices satisfying the strong 
ellipticity condition, i.e., there is a number A > such that 

N n 

(2.2) <(^e>A||| 2 :=A]T£|a| 2 , VxGM" 

i—1 a— 1 
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We also assume that Aff are bounded, i.e., there is a number A > such that 

N n 

(2.3) Kf (*)l 2 < A 2 , Vxet" 

i,j=l q,/3— 1 



If we write (|2.1j) component-wise, then we have 

(2.4) (Lui^-D^AgDfiu 1 ), Vi = l,...,N. 
The transpose operator of 'L of L is defined by 

(2.5) l Lu = -D a ?A ap D fi u), 

where 'A Q/3 = (A /3q ) t (i.e., 'A?f = A??). Note that the coefficients *A a / satisfy 
(12.21). (12.311 with the same constants A, A. 



In the sequel, we shall use the notation f s f :— J s f (assuming < |5| < oo), 
where S is a measurable subset of R™ and |5| denotes the Lebesgue measure of 
measurable S. 

Definition 2.1. We say that the operator L satisfies the property (H) if there exist 
Ho, Hq > such that all weak solutions u of Lu = in Br = Br(xq) satisfy 



2 /T\ n ~ 



B r 



(2.6) / \Du\'<H [-) / \ Du \ 2 i 0<r<s<R. 



Similarly, we say that the transpose operator *L satisfies the property (H) if corre- 
sponding estimates hold for all weak solutions u of l Lu = in Br. 

Lemma 2.2. Let (a Q/3 (a;))™ p =1 be coefficients satisfying the following conditions: 
There are constants Ao , Ao > such that for all x £ M n 

(2.7) a afi (x)i^ a > Ao|C| 2 , V£eK n ; ]T \a aP (x)\ 2 < A 2 . 
Then, there exists eo = eo( n , Ao, Ao) such that if 

(2.8) e 2 {x) := £ |A<f (z) - a Q/3 (:r)<% |* < e 2 , g R", 

t/ien £/ie operator L associated with the coefficients A^ satisfies the condition (H) 
with Ho = Ho(n, A , A ),-ff = #o(n, AT, A , A ) > 0. 

Proof. See e.g., PU Proposition 2.1]. □ 

Lemma 2.3. Suppose that the operator L satisfies the following Holder property 
for weak solutions: There are constants Ho, Co > such that all weak solutions u 
of Lu = in Bir = B-2r{xq) satisfy the estimate 



(2.9) 



[u] c » {b r ) <C R-»° |u 



1/2 

2 



where [f]c(Q) denotes the usual C*(f2) semi-norm of f : see [lOj for the definition. 
Then, the operator L satisfies the property (H) with ho & n d Hq = Ho(n, N, A, A, Co). 
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Proof. We may assume that r < s/4; otherwise, (|2.6|) is trivial. Denote u r — f B u. 
We may assume, by replacing u by u — u s , if necessary, that u s = 0. From the 
Caccioppoli inequality, (|2.9p . and then the Poincare inequality, it follows 

\u(x) — u(y)\ 2 dydx 



f \Du\ 2 <Cr- 2 f \u-u 2r \ 2 <Cr~ 2 f j 



< Cr- 2 [u} 2 cm{B2r) (2r) 2 ^ \B 2r \ < Cr n - 2+2 »°[u] 2 CM{Bs/2 

<C{r/s) n - 2+2 ^s- 2 f \u\ 2 <C(r/s) n - 2+2 ^ f \Du 
Jb s Jb b 



2 



(2-12) \W\ L ~(B R/2 )<C\f M 



The proof is complete. □ 

Lemma 2.4. Assume that the operator L satisfies the property (H). Then, the 
operator L satisfies the Holder property (|2.9[) . Moreover, for any p > 0, there 
exists C p = C p (n, N, X,A, no, Hq,p) > such that all weak solutions u of Lu = 
in Br — Br(xq) satisfy 

(2-10) ll«llz- (Br ) < {R Cp r)n/p Mi*(B») > Vr e (0, R). 

Proof. From a theorem of Morrey (TTJ, Thoerem 3.5.2], the property (H), and the 
Caccioppoli inequality, it follows that 

(2.11) [u] 2 cn(BR) < C R 2 - n - 2 ^ \\Du\\ 2 L2(B3R/2) < CR- n - 2 ^ \\u\\ 2 l2(B2r) . 
Then, by a well known averaging argument (see e.g., [12] ) we derive 

1/2 

2 ' 

\U 

>B 2 R 

where C = C(n, N, A, A, /i , H ) > 0. For the proof that (f2~T2^) implies f2~TD]) . we 
refer to pp. 80-82]. □ 

2.2. Function spaces Y 1 > 2 (Q) and Y^' 2 ^). 

Definition 2.5. For an open set SI C R n (n > 3), the space Y 1,2 (ti) is defined as 
the family of all weakly differentiable functions u E L 2 (O), where 2* = -^g, whose 
weak derivatives are functions in L 2 (fl). The space Y 1,2 (Sl) is endowed with the 
norm 

IM!yi.2(S2) : = ll u IL=*(n) + \\Du\\ L2{n) . 

We define F 1,2 (O) as the closure of C™(fl) in Y^ 2 (Q), where C£°(0) is the set of 
all infinitely differentiable functions with compact supports in SI. 

We note that in the case SI = R", it is well known that Y 1 ' 2 (R n ) = F 1,2 (R n ); 
see e.g., [El p. 46]. By the Sobolev inequality, it follows that 

(2-13) \\u\\ L2 * (n) <C(n)\\Du\\ L2{n) , Vu E Y^ 2 

Therefore, we have Wq' 2 (SI) C Y h2 (Sl) and Wq' 2 (SI) = Y h2 (Sl) when SI has finite 
Lebesgue measure. 

From (|2.13|) . it follows that the bilinear form 

(2.14) (u,v) u := / D a u l D a v l 
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defines an inner product on H := Y ' (Q) N . Also, it is routine to check that H 
equipped with the inner product (12.14|) is a Hilbert space. 

Definition 2.6. We shall denote by H the Hilbert space YQ' 2 (n) N with the inner 
product (|2.14p . We denote 

IMIh : = ( w > m )h 2 = ll-^HIz^n) • 
We also define the bilinear form associated to the operator L as 



B{u,v) := / AfDpuW^ 
Jn 



By the strong ellipticity (|2.2[) . it follows that the bilinear form B is coercive; i.e, 
(2.15) B{u,u) > X(u,u) n . 

3. Fundamental matrix in R" 

Throughout this section, we assume that the operators L and *L satisfy the 
property (H). The main goal of this section is to construct the fundamental matrix 
of the the operator L in the entire R™, where n > 3. Since F 1,2 (R") = F 1,2 (R"), 
we have, as in Definition I2J51 



llull^-CR-) < C(n) ||£>u|| ia(RB) , Vu G F^(R"). 

We note that W 1 ' 2 (R n ) C Y 1 ' 2 (W l ) C W^(W l ). Unless otherwise stated, we 
employ the letter C to denote a constant depending on n, N, A, A, p Q , H Q , and 
sometimes on an exponent p characterizing Lebesgue classes. It should be under- 
stood that C may vary from line to line. 

3.1. Averaged fundamental matrix. Our approach here is based on that in 
[IT] . Let y G R" and 1 < k < N be fixed. For p > 0, consider the linear functional 
f D , sU k . Since 

J B P (y) 



u k 



(3.1) 

<B p (y) 

Lax-Milgram lemma implies that there exists a unique v p = v p . y ^k G H such that 
(3.2) f Al P D p vlD a u l = / u k , VueH. 

Note that (f27T5|) . |3?2| . and ([3~Tj) imply that 

A||^||^<B(« P ,^)<Cp( 2 -")/ 2 ||t; p || H , 

and thus we have 

(3-3) \\Dv p \\ L2(Mn) = \\v p \\ H <Cp^ 2 . 

We define the "averaged fundamental matrix" T p ( ■ ,y) = (r? fc ( • , y))f k=1 by 

(3-4) I* (.,») = «2=«J. yifc . 

Note that we have 

(3.5) / A^D f3 T P k {-,y)D a u^i u k , Vtt G H, 

A" ./l? p (3/) 
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and equivalently (a <-> f3 , i j). 

(3.6) / t AfD uW a T p ik (-,y) = -f u\ Vit G H. 

JB» "'S p (y) 

In the sequel, we shall denote by L£°(f2) the family of all L°° functions with 
compact supports in O. For a given / £ L^R™)^ consider a linear functional 

(3.7) w h-> / / ■ w, 

JR™ 

which is bounded on H since 



(3.8) 



< 11/11 _2a_ llwll Jt <C 11/11 _3n_ \\w\\„ 



Therefore, by Lax-Milgram lemma, there exists m g H such that 

(3.9) / t A&D ft v?D a vf = [ fw\ VweH. 

JR" JR" 

In particular, if we set w = v p in (|3.9p . then by (|3.6p . we have 

(3.10) / re k (.,y)f = -f u k . 

Jr™ Jb p ( v ) 
Moreover, by setting w = u in (|3.9p . it follows from (|3.8p that 
(3.H) ||-D U || i2(Rn) <C||/l| i2n/(n+2)(Rn) . 

3.2. L°° estimates for averaged fundamental matrix. Let it G H be given as 
in (|3.9p . We will obtain local L°° estimates for u in Br(xq), where xq £ M. n and 
R > are fixed but arbitrary. 

Fix x £ Br(xo) and < s < R. We decompose u as u — U\ + U2, where 
ui £ W 1 ' 2 (B s (x)) N is the weak solution of l Lu\ = in B s (x) satisfying u\ — u on 
dB s (x); i.e., Ui — u £ W ' (B s (x)). Then, for < r < s, we have 

\Du\ 2 < 2 f \Du x \ 2 + 2 j \Du 2 \ 2 

|2 



B r {x) JB r (x) JB T (x 



<C(- / \Dm\ 2 + 2 I \Du 2 

" s/ Jb s {x) Jb 3 (x) 

<C[-) / \Du\ 2 + C \Du 2 \ 2 . 

Vs/ Jb 3 (x) Jb s (x) 

Since u 2 £ Wq' 2 (B s (x)) n is a weak solution of t Lu 2 = / in B s (x), we have 
/ \Du 2 \ 2 <C\\f\\ 2 L2n/(n+2){B {x)) . 

JB s (x) 

For given p > n/2, choose po £ (n/2,p) such that fj,x := 2 — n/po < Mo- Then 
(3.12) 11/11^^ < ll/llk(B s{!C ))l^| 1+2/ ^ 2/P0 <C||/||L 0(an)S »- 2+2 «. 
Therefore, after combining the above inequalities, we have for all r < s < R 



f f \ n — 2+2un 

|D«r <C(- 



/ |^tt| 2 + Cs"- 2+2w ||/| 

JBJx') 



I 2 

Ilpo(R«) 

'B r (sc) Ni/ JB s (x) 
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By a well known iteration argument (see e.g., Lemma 2.1, p. 86]), we have 
/ \Du\ 2 <C{-) / |£«| 2 + Cr- 2 + 2 ^ H/lli™^) 

/ r s n-2+2^i C „ 

< c (-) y Rn |du| 2 + cv- 2 + 2 ^ n/ii^ (R , l} , 

for all < r < R and x £ B R (x a ). From (13.13j) it follows (see, e.g. [H]) 

(3-14) Ncmx(B r( xo)) < c(i?-("- 2+2 ^) IpuH^^) + H/llLoCRn)) • 

Note that since it 6 H, we have 

IMI^so)) ^ II^II'-cbh^)) l^l 2/ " < ci? 2 pti||* a(RB) . 

Consequently, we have 

\\ U Wl°°(B r/2 (x )) ^ ^^'^C^B^o)) + C '^ " \\ U \\l2(B r (x )) 

< C (i? 2 "" \\Du\\ 2 L2{Mn) + R 2 ^ ||/||Lo(Rn)) + CR 2 - n \\Du\\ 2 L2(M „ 



< CR 2 " ||/||£2n/(n+2)(R») + CR 2 ^ 1 ||/||ip (Rn) l 

where we used the inequality (|3.11[) in the last step. 

Therefore, if / is supported in Br(xo), then (|3 . 1 2[) yields (recall fi\ = 2 — n/po) 

(3-15) \\u\\ l ^ Br/2{xo)) < CR 2 -'^ \\f\\ LPO{Ba{xo)) < CR 2 -^ \\f\\ LP{BR{xo)) • 

Now, (|3.10p implies that for p < R/2, we have, by setting x n = y in (|3.15p . 



(3.16) 



</ I«I<^ 2 ^ /p II/IIl, (Sj? (,)), Vp>n/2 

" , -B p (j / ) 



provided that / is supported in B R (y). Therefore, by duality, we see that 

(3.17) \\v p \\ Lq{BR[y)) <CR 2 - n+n l\ Vge[l,^), Vpe(0,fl/2), 

where v p = v p . y ^ is as in (|3.4| . 

Fix x ^ y and let r := | |ar — y|. If p < r/2, then since v p £ W 1 ' 2 (B r (x)) N and 
satisfies Lv p = weakly in B r (x), it follows from Lemma \2. 41 that 

(3.18) \v p (x)\ < Cr- \\v p \\ L1{BAx)) < Cr- \\v p \\ L1(Bar(y)) < Cr 2 - . 
Since p, y, k are arbitrary, we have obtained the following estimates. 

(3.19) \T"(x,y)\ <C\x-y\ 2 - n , Vp < \x - y\ /3. 

3.3. Uniform weak-L^ 5 estimates for T p ( - ,y). We claim that the following 
estimate holds: 



(3.20) / |r p ( •, y)\— < CR~ n , VR>0, Vp > 0. 

JR™\B R (y) 

HR> 3p, then by (f3~19|) we have 

f \T p {x,y)\^ dx <C [ \x-y\~ 2n dx <CR- 

«"\Bh(9) JWl™\B R ( y ) 
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Next, we consider the case R < 3p. Let v T p be the fc-th column of the averaged 
fundamental matrix T p (- ,y) as in (|3 .4|) . From (|3.3[) . we see that 

\\ v p\\l 2 ' (R™\B R {y)) - W V p\\l 2 * (&") - W Dv p\\l 2 (R") - ^ P^ " )/2 - 

and thus (|3.20[) also follows in the case when R < 3p. 

Now, let A t = {x £ R™ : |r p (ir, y)\ > t] and choose R = £-V("-2). Then, 

\A t \B R (y)\<t-^ [ \T"(.,y)\^ <Ct-^t^ =Ct~^. 

JA t \B R {y) 

Obviously, \A t <~]B R (y)\ < CR n = Ct~^>. Therefore, we obtained that for all 
t > 0, we have 

(3.21) \{x S R" : \T p {x lV )\ > t}\ < Gt~£* , Vp > 0. 

3.4. Uniform weak-i"^ 1 estimates for DT p (-y). Let v p be as before. Fix a 
cut-off function rj 6 C°°(R™) such that r\ = on B R / 2 (y), r) = 1 outside B R (y), 
and ID77I < C/i?. If we set it := ry 2 f (9 , then by 



= / V 2 AfD v p D a v p + I 2 V AfD^ p v;D aV , 
which together with (|3.19[) implies that if R > 6p, then 



\Dv P Y < CRT \ \v p \ z < CR 2 

' J flsW\B«/a{v) 

On the other hand, if i? < 6p, then (|3 . 3|) again implies 

/ |£^l 2 < / I^pI 2 < c P 2 ~ n < CR 2 - n . 

JM n \B R (y) ' JR n 

Therefore, we have 

(3.22) f \DT p (- iy )\ 2 < CR 2 - n , Vi? >0, Vp > 0. 

JR™\B R (y) 

Next, let A t = {x e R" : \D x T p (x, y)\ > t} and choose R = t-V(n-i). Then 
\At\B R (y)\<r 2 [ \DT p (-,y)\ 2 <Ct-^ 

JA t \B R (y) 

and |A t ("1 < CR n = Ct~^. We have thus find that for all t > 0, we have 

(3.23) I {a; 6 R n : \D x T p (x, y)\ > t}\ < Ct~^, Vp > 0. 

3.5. Construction of the fundamental matrix. First, we claim 

(3.24) \\DT p {-,y)\\ LP[BR(y)) <C p R l - n+n l p , Vp > 0, Vp G (0, ^). 
Let v p be as before. Note that 

\Dv p \ p = [ \Dv p \ p + [ \Dv p \ p 

B R (y) J B R (y)n{\DvP\<T} J B R (y)n{\DvP\>r} 

<t p \B r \+ [ \Dv p \ p . 

J{\DvP\>t} 
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By using ()3.23|) . we estimate 

f \Dv p \ p = f pt^ 1 \{\Dv p \ > max(i,r)}| dt 

J{\Dvp\>t} Jo 



< Ct 



pt p - l dt + C I piP-W(»-i) dt 



cfl-p/ip-^rr-^-V. 



By optimizing over r, we get 
(3.25) 



/ \Dv p \ p < CR {1 ~ n)p+n , 



,v) 



from which (|3.24[) follows. 

If we utilize (|3.21[) instead of (|3.23[) . we obtain a similar estimates for T p ( 

(3.26) \\r p (-,y)\\ LHBR(y)) <C p R 2 - n+n ^, V P >0, Vpe(0,^). 

Let us fix q € (1, — zr)- We have seen that for all R > 0, there exists some 
C(R) < oo such that 

\\r p (-,y)\\w^(Bn(y)) < C W. v / 5>0 - 

Therefore, by a diagonalization process, we obtain a sequence and r( • , y) 

in W r ; 1 oc 9 (R n ) AfxA ' such that lim^oo = and that 



(3.27) 



r""(.,y)^r(.,») in W^(B R (y)) 



NxN 



Vi? > 0, 



where we recall that — 1 denotes weak convergence. Then, for any cj) € (M. n ) N , 
it follows from (13.51) 



(3.28) 



/ AfD p T jk ( • , y)D a tf = hm / A^D p Y% ( • , y)D a < 

JR™ fi^>oo J Rn 



lim f r Afe 



\y)- 



Let v T p be the fc-th column of T p ( ■ , y) as before, and let i> T be the corresponding 
jfe-th column of T( • , y). Then, for any g G Lf(B R (y)) N , yields 



(3.29) 



/ « 9 


= lim 


/ • 




/i — >oo 


JR™ 



< C p i? 2 -™+"/Hl9|| LP ' (Bfl(!/)) 



where p' denotes the conjugate exponent of p € [1, z^zn)- Therefore, we obtain 

(3.30) \\T{-,y)\\ LHBR{ y))<C P R 2 - n+n/p , V P e[i,^). 

By a similar reasoning, we also have by (|3.24[) 

(3.31) \\DT{-,y)\\ LP[BR{y)) <C p R}- n+n lr, \/ p€ [l,-^). 
Also, with the aid of ([3T20"]) and (|3?22| . we obtain 



(3.32) 
(3.33) 



n \B R (y) 



\T(-,y)\ 2 <CR~ n , 



l \B R (y) 



\DT(-,y)\ 2 <CR 



2-n 
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In particular, ()3.32|) . (|3.33|) imply that 

(3-34) l|r(-,y)|| rl , 2(RnSr(t/)) <Cr 1 -"/2 ! Vr>0 . 

Moreover, arguing as before, we see that the estimates (13.32[) and (|3.33[) imply 

(3.35) : \T(x,y)\ > t}\ <Ct~^, V<> 

(3.36) \{x e R™ : \D x T(x,y)\ > t}\ < Ct~^ W > 0. 



Next, we turn to pointwise bounds for r( • , y). Let v T be the fc-th column of 
r( • , y). For each x ^ y, denote r — | \x — y\. Then, it follows from (|3.34|) and 
(I3.28P that v is a weak solution of Lv = in B r (x). Therefore, by Lemma |2"^31 and 
(ESP we find 

(3.37) \v(x)\ < Cr~ n ||«|| £l(Bp(a)) < O"" \\v\\ L1{Bar{y)) < Cr 2 ~ n , 
from which it follows 

(3.38) \T(x,y)\<C\x-y\ 2 - n , Vx^y. 

3.6. Continuity of the fundamental matrix. From the property (H), it follows 
that r( • , y) is Holder continuous in R" \ {y}. In fact, (|2~TTj) together with (f3~2"g)) 
and (|3.33j) implies 

(3.39) \T(x,y)-T(z,y)\<C\x-z\ t *°\x-y\ 2 - n -> 10 if \x - z\ < \x - y\ /2. 

Moreover, by the same reasoning, it follows from (|2.1ip and (|3 . 22[) that for any 
given compact set K <£ 1" \ {y}, the sequence {F Pl * ( ■ , y)}°^ =1 is equicontinuous on 
K. Also, by Lemma [2.41 and (|3.20|) . we find that there are Ck < oo and pk > 
such that 

(3.40) ||r"( • , y)\\ LBe(K) < C K Vp < p K for any compact K € K" \ {y} . 
Therefore, we may assume, by passing if necessary to a subsequence, that 

(3.41) T Pfi ( ■ , y) — > T( • , y) uniformly on K, for any compact K <s R™ \ {y} . 

We will now show that T(x, ■ ) is also Holder continuous in R™ \ {x}. Denote 
by T" T ( ■ , x) the averaged fundamental matrix associated to t L, the transpose of L. 
Since each column of T p ( ■ , y) and T" T ( ■ , x) belongs to H, we have by (|3.5[) . 



(3.42) 




rf fc (-,y)- 



By the same argument as appears in Sec. 13.51 we obtain a sequence {<J^}^ =1 tending 
to such that i T a " ( • , x) converges to *T( • , x) uniformly on any compact subset of 
W 1 \ {x}, where *T( ■ ,x) is a fundamental matrix for *L satisfying all properties 
stated in Sec. S3] By ([£12"]) . we find that 

nr(->*) = / r «(-.y)- 

From the continuity of T^ 1 ( • , y), it follows that for x, y G 1" with x ^ y, we have 
lim ff « = lim / r^(-,y)=rf fe -( a; ,y) 
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and thus by (|3.41j) we obtain 

lim lim = lim T^(x,y) = T lk (x,y). 

On the other hand, (|3.27|) yields 

lim ff « = lim / <T%(-,x)=-f 

and thus it follows from the continuity of *Tm( • , x) that 

lim lim 5" = lim / t T k i(- ,x) = t T k l(y,x). 

We have thus shown that 

Ti k (x,y) = *r fc ;(y,x), Vfc, I = 1, . . . , N, Vx^y, 
which is equivalent to say 

(3.43) T(x,y)^ t T(y,x) T , Vx^ty. 

Therefore, we have proved the claim that T(x, ■ ) is Holder continuous in R™ \ {x}. 

So far, we have seen that there is a sequence {p^Y^—i tending to such that 
rPM ( • , V) -> r( • , y) in R" \ {y}. However, by ([3~42]) . we obtain 

Tf fe (x, y) = lim / rf fe ( ■ , y) = lim / T% ( ■ , x) 

(3 44) ^Jb^x) U ^°°JB p (y) 

= / t r kl (-,x) = -f r lk (x,-), 

JB p (y) JB p (y) 

i.e., we have the following representation for the averaged fundamental matrix: 

(3.45) T p (x,y) = f T(x,z)dz. 

JB p {y) 

Therefore, by the continuity, we obtain 

(3.46) ]imT"{x,y) = T{x,y), x^y. 

3.7. Properties of fundamental matrix. We record what we obtained so far in 
the following theorem: 

Theorem 3.1. Assume that operators L and L satisfy the property (H). Then, 
there exists a unique fundamental matrix T(x, y) = (Tij(x, y))fj =1 (a; ^ y) which is 
continuous in {(x,y) G 1" x R™ : x ^ y} and such that T(x, ■ ) is locally integrable 
in 1" for all x S R" and that for all f = (/\ . . . , f N ) T G C™(R n ) N , the function 
u = (u 1 , . . . , u N ) T given by 

(3-47) u(x) := f T(x,y)f(y)dy 

belongs to Y 1,2 (R n ) N and satisfies Lu — f in the sense 

(3.48) J AffD u j D a ^ = J f<j>\ V0eC c oo (R") A '. 
Moreover, T(x,y) has the property 

(3.49) / AffD p T jk (-,y)D a ^ = k (y), V0 e C™(R n ) N . 
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Furthermore, T(x,y) satisfies the following estimates: 

(3.50) ||r< • . v)\\y^\ba v )) + H r (^ ' )llr^(K»\fl.(x)) < Cr 1 -^ , Vr > 0, 

(3.51) ||r( • , y)\\ LP{Briy)) + \\T(x, ■ )\\ LP{Br{x)) < C p r 2 - n+ ^ , Vp G [1, 

(3.52) \\Dr(-,y)\\ LP{BAy)) + \\DT(x, ■ )\\ LP(Br{x)) < C^r 1 ""*?, Vp G [1, 

(3.53) |{z G 1" : |r(a;,y)| > t}\ + \{y G M" : \T(x,y)\ > t}\ < Ct~&, 

(3.54) \{x G R" : |DJ(x,i/)| > t}| + \{y G R" : |-D„r(;r, y)| > t}| < Cf^, 

(3.55) |r(x,y)| <C|x-y| 2 -", Vx ^ y, 

(3.56) |r(x,y)-r(z,y)| <C|x-zr Ix-yl 2 -"^ |ar - z| < |x - y\ /2, 

(3.57) ir^yJ-Tfx^)! <C|y-^r \x-y\ 2 - n ~^ if \y - z\ < \x - y\ /2, 

where C — C(n, AT, A, A, ^,q,Hq) > and C p = C p (n, N, A, A, /xo, Hq,p) > 0. 

Proof. Let r p (x,y) and r(x, y) be constructed as above. We have already seen 
that r is continuous in {(x, j/) G 1™ x R™ : x ^ y} and satisfies all the properties 
(|3.49[) - (|3.57[) . By using the Lax-Milgram lemma as in Sec. 13. li we find that for 
all / G C™(R n ) N , there is a unique u G F 1 ' 2 (R n ) w satisfying 

f AffDpuWaV^ f /V, Vv G Y 1,2 (R n ) w . 

If we set i/ = r?.(x, • ) above, then (|3 . 5|) together with (|3.43[) implies that 

(3-58) / T^(x,-)f=[ t A^D a ^ ik (.,x)D^ = -f u k . 

Assume that / is supported in Br(x) for some R > 0. Then, by (|3.27[) and (|3.43|) 
we have 

lim / T^ix, ■ )f = lim / T p ki (x, ■ )f = f T kl (x, ■ )f\ 

By the same argument which lead to (|3.14[) in Section |3~!?1 we find that u is Holder 
continuous. Therefore, (|3.47[) follows by taking the limits in (|3.58p . 

Now, it only remains to prove the uniqueness. Assume that T(x, y) is another ma- 
trix such that r is continuous on {(a;, y) G R™ x R™ : x ^ y\ and such that T(x, ■ ) 
is locally integrable in R™ for all x G R™ and that for all / G C%°(M. n ) N , 

u(x) := / t(x,y)f(y)dy 

belongs to Y" 1,2 (IR™) and satisfies Lu = f in the sense of (I3.48p . Then by the 
uniqueness in H = Y 1 ' 2 (M n ) N , we must have u = u. Therefore, for all x G R™ we 
have 

/ (T-F)(x,-)f = 0, V/6Cr(KT, 
and thus we have T = f in {(x, y) £ R" xl":i / y}. □ 

Theorem 3.2. Assume that the operators L and L satisfy the property (H). If 
f G (L >»+ 2 (R") n Lf (W n )) N for some p > njl, then there exists a unique u in 
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Y 1 ' 2 (R n ) N such that 

(3.59) f AffD uW a v l = [ /V, Vw 6 y^QET)* 
Moreover, u is continuous and has the following representation: 

(3.60) u k (x)= [ T ki (x,y)f(y)dy, k=l,...,N, 



where (T k i(x,y)) k i=1 is the fundamental matrix of L. 

Proof. Since / G L™^(W. n ) N , the same argument as appears in Sec. 13.11 implies 
that there is u e F 1 » 2 (R n ) JV satisfying ([3"39")> . If we set v 1 = T p ki (x, •) in 1(539]) . 
then (|3.2[) implies that 

(3-61) / r p ki (x,.)f= [ t A^D a t T p ik (-,x)D^=-f „ 

Next, note that ((3~20l) . (029, and the assumption f £ L^(R n ) n £f oc (R n ) for 
some p > n/2, imply that 

lim / i* (*, ■ )f = iim [ / r*(*, ■ )f + [ r p ki (x, ■ )A 



A 



Bi(x) JH»\Bi(x) 



( 3 - 62 ) = / r fcl (x, •)/*+/ r fci (z, .)f 

JB 1 (x) Jw. n \B 1 (x) 

r ki (x, -)f. 

Finally, by the same argument which lead to (|3.14p in Sec. 13. 2[ we find that u is 
Holder continuous, and thus (|3.60|) follows from (13.61)1 and (I3.62[) . □ 

Corollary 3.3. Suppose that f — (f 1 , . . . , f N ) T has a bound 

(3.63) |/(x)| <C(l + \x\)- (1+n/2+e) Va; € R" for some e > 0. 

Then, u — (it 1 , . . . , u N ) T given by p.60[) is a unique Y 1 ' 2 (W n ) N solution of Lu = f 
in R™ in the sense of (I3.59p . 

Proof. Note that (|3l%|l implies / E (i^(R n ) n L p (R")) iV '. □ 

Theorem 3.4. Assume that L and *L satisfy the property (H). If f 6 F 1,2 (R") JV 
satisfies Df G -Lj' (R™)^*™ /or some p > n, then 



(3.64) f k (x) = J D a T kl (x, ■ )A% D [3 p, k=l, 



,N, 



where (T k i(x, y)) k i=1 is the fundamental matrix of L. 

Proof. We denote by *T P the averaged fundamental matrix of t L. Recall that 
columns of *T P belong to H. Then, by (|3 . 5[) we have 



/ t A^D a t T P k (-,x)D p f = i f k 

JR" JB p {x) 
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As in ([3J2]) . the assumption Df G Lf oc (R n ) N for p > n, together with (f3~22]) and 
(f3~24j) yields 

lim / f = lim ( / + / 'AfrDcXP,, ( • , x)D p A 

JBi(x) JR"\B 1 (x) 



f AfD a T kl (x,-)Dpf 



where we used (|3.43p in the last step. By the Morrey's inequality [17], / is contin- 
uous and thus (|3.64|) follows from (|3.65|) . □ 



Corollary 3.5. Assume that L, l L, L, and L satisfy the property (H). Denote by 

,a/3 



r and r the fundamental matrices of L and L, respectively. If the coefficients A" 



of L and Afj 3 of L are Holder continuous, then 

(3.66) t lm (x,y) = T lm (x,y) + / D a T h (x, ■ )(A°f - A^ )Dpf jm ( ■ , y), x^y. 

Proof. We denote by T p and T p (p < \x — y\ /4) the averaged fundamental matrices 
of L and L respectively. Recall that columns of T p and T p belong to H. Moreover, 
since we assume that the coefficients are Holder continuous, the standard elliptic 
theory, (j3~38"]) . and implies that DT p (x, ■ ) and DT P ( ■ , y) are locally bounded. 

Therefore, by setting / J = T p m ( ■ , y) in (|3.64[) we have 

(3.67) ff m 0r, y) = / D a T H (x, ■ )A«f D f p m (- , y), 
Next, set / J = F P j(x, ■ ) and apply (|3.64[) with L replaced by *L to get 

rf„>, V)= I DjT^y, ■ yAfDpT^x, ■ ). 

JE" 

By using (|3.43p and interchanging indices (a /3, i j), we obtain 

(3.68) T p m (x, y) = [ D a T» u {x, ■ )A%?Dpf jm {- , y). 



Now, set r = \x — y\ /A and split the integral (|3.67|) into three pieces (recall p < r) 



D a T H (x,.)A^ D T» jm {.,y). 

B r (x) JB,.(y) JR n \(B r (x)llB r (x)) 

Since we assume that the coefficients are Holder continuous, it follows from the stan- 
dard elliptic theory that DT(x, ■ ) and DT( ■ , y) are continuous (and thus bounded) 
on B r (y) and B r (x) respectively. Moreover, (|3.45[) implies 

Df p ( • ,y) -> DT( ■ , y) uniformly on B r (x) as p -> 0. 

Therefore, as in (|3.65[) , we may take the limit p — > in (|3 .67[) to get 

Tim{x, y)= j D a T h (x, ■ )A c *?Dpf jm (- , y), 
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Similarly, by taking the limit p — ► in (|3.68| . we obtain 

Fi m {x,y)= / D a Tu(x, ')AffDpfj m (- ,y). 



3 

The proof is complete. □ 
Remark 3.6. We note that in terms of matrix multiplication (|3.60p is written as 

u(x) = / T(x,y)f(y)dy, 

JR™ 

where both u, f are understood as column vectors. Also, (|3.66p reads 
f (s, y) = T(x, y)+ [ D a T( Xl ■ ){A a ? - A a ?)DpT{ ■ , y). 

JR™ 

4. Green's matrix in general domains 

4.1. Construction of Green's matrix. In this section, we shall construct the 
Green's matrix in any open, connected set C R™, where n > 3. To construct the 
Green's matrix in fl, we need to adjust arguments in Section [3] 

Henceforth, we shall denote £l r (y) '■= ^ H B r (y) and d y := dist(j/, dVi). Also, as 
in Section[3l we use the letter C to denote a constant depending on n, N, A, A, /xo, 
Hq, and sometimes on an exponent p characterizing Lebesgue classes. 

It is routine to check that for any given y G fi and 1 < k < N, the linear 
functional u i— > f n . , w fc is bounded on H = ig ' (Q)^. Therefore, by Lax-Milgram 
lemma, there exists a unique v p = v pv y t k G H such that 

(4.1) / AfD v p D a u l = S u\ VttGH. 

Note that as in (13.31). we have 



(4-2) \\Dv p \\ L2m = \\v p \\n<C\£l p {y)\-*r. 

We define the "averaged Green's matrix" G p ( ■ , y) = (G P k ( ■ , y))f k=1 by 

G p jk {-,y) = v 3 p = v p . Vtk . 

Note that as in (|3 . 5[> . we have 

(4.3) [ A^D G p (-,y)D a u l = -f u\ Vtt G H. 

Next, observe that as in (|3.7j) — fj3. 10|) . for any given / G L^°(Q) N , there exists a 
unique u G H such that 

Pi. «A f « 1 u fe 



G? fc (-,y)f = 

Moreover, as in (|3 . llj) . we have 

\\Du\\ L 2(ty < C ||/||i2n/(n+2)(n) • 

Also, by following the argument as appears in Section 13.21 we find that if / is 
supported in Bn(y), then we have 

H\^ (BBMy)) < CR 2 ^ \\f\\ LP(BR(y)) , Vi? < d y , Vp > n/2. 
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Therefore, as in (|3.16|) . for any / £ L£°(Br(2/)), R < d y , we have 

G p lk (-,y)f <CR 2 - n /P\\f\\ LP{Bniy)V Vp<i?/4, Vp>n/2. 

Bn(y) 

Therefore, as in (|3.17p . we see that if R < d y , then 

\\G p {-,v)\\ Lq{BR(y)) <CR 2 - n+n l\ Vp<i?/4, Vqg[1,^2). 
Then, by following the lines in (|3.18[) - (|3.19p . we obtain 



\G p (x,y)\ <C\x-y\ 2 ~ n if \x - y\ < d y /2, Vp<\x-y\/3. 

Next, we shall derive an estimate corresponding to (|3.22|) . Let 77 £ C°°(K n ) be 
a cut-off function such that < 77 < 1, r/ = l outside B R / 2 (y), T] = on B R /4(y), 
and \Drj\ < C / R, where R < d y . By setting u = rj 2 v p £ H in (|4.ip . we obtain 



f v 2 \Dv p \ 2 <C f \D V \ 2 \v p \ 2 < CR- 2 f 
Jn Jn Jb„ 



|2 
\Vp\ 

/2(y)\B R/i (y) 



( 4 - 4 ) <CR- 2 [ \x-y\ 2 ^dx 

J Bn/2(y)\B R /4(y) 

= CR- 2 R 4 - n = CR 2 -", Vp < R/12. 
Therefore, we have (r = R/2) 

(4.5) f \DG P {- ,y)\ 2 <Cr 2 - n 7 Vp < r/6, Vr < dy/2. 
On the other hand, (14. 2p implies that if p > r/6, then 

(4.6) / \DG"(-,y)\ 2 < [ \DG?{-,y)\ 2 <C\tl p {y)\^ <Cr 2 -. 
Jn\B r (y) Jn 

Therefore, by combining (|4. 5[) and (|4.6p . we obtain 

(4.7) / \DG p (-,y)\ 2 < Cr 2 - n , Vr < dy/2, Vp > 0. 
Jn\B r (y) 

From the estimate (|4.7I) , which corresponds to (|3.22p , we can derive an estimate 
corresponding to (|3.24p as follows. By following the lines between (|3.22|) and (|3.23p . 
we obtain 

(4.8) \{x£Q.:\D x G p (x,y)\>t}\<Cr^, Vp > if t > (d y /2) 1 - n . 
Then, by following lines (|3~34) - (l3~2"51) . we find (set r = (i?/2) 1 -™) 

(4.9) / \DG"(-,y)\ p <CW^ +n , VR < d y , Vp > 0, Vp G (0, 

Now, we will derive estimates corresponding (|3.20p and (|3.26|) . Let 77 be the 
same as in (|4.4p . Note that (14. 4p and (|4.7p implies that for R < d y , 

(4.10) / |i?(^p)| 2 < 2 / /7 2 |£v p | 2 + 2 / \Dr)f \v p \ 2 < CR 2 - n , Vp < R/12. 
Jn Jn Jn 

Since r/v p £ H = y Q 1:2 (^), it follows from (jCTJ]) and (|2~ilij) that 

(4.11) / \v p \ r <Cr~ n , Mr < d y /2, Vp < r/6. 

jQ,\B r {y) 
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On the other hand, if p > r/6, then (14. 2| implies 



(4.12) 




<C\9. p \~ l < Cr~ 



Therefore, by combining (|4.11[) and (|4.12[) . we obtain 

(4.13) f \G p (-,y)\ 2 ' <Cr~ n , Wr<d y /2, Vp > 0. 

Jn\B r (y) 

As in Section l3~3l the above estimate (|4.13|) yields 

(4.14) \{x e fi : \G"{x,y)\ > t}\ < Ct~^ , Vp > if t > {d y /2) 2 ' n . 
Then, as we argued in (|4.9p . we find (set r = (i?/2) 2 ~") 

(4.15) / \G»(-,y)\ p <CBP^+ n , VR<d y , Vp > 0, Vp G (0, ^). 

Now, observe that (|4. 9|) and (|4.15p in particular imply that 

(4.16) \\G P ( - ,y)\\ w i, P{Bdy(y)) < for somep G (1,^), uniformly in p. 

Therefore, from (|4.16|) together with (|4. Tjl and (|4.13p . it follows that there exist a 
sequence {p^Y^i tending to and functions G( ■ , y) and G( ■ , y) such that 

(4.17) y)^G(-,y) in W^(B dy (y)) N * N and 

(4.18) G ?»(-,y)^G(-,y) in F 1 * 2 ^ \ B dy/2 (y)) N * N as p -> 00. 

Since G( ■ ,y) = G( ■ , y) on (y) \ B dy/2 {y), we shall extend G( • , y) to entire ft 

by setting G( • , y) = G( ■ , y) on fl \ Bd y (y) but still call it G( • , y) in the sequel. 
Moreover, by applying a diagonalization process and passing to a subsequence, if 
necessary, we may assume that 

(4.19) G p »{-,y)^G{-,y) in y J ' 2 (fi \ B r (y)) NxN as /1 -» 00, Vr < 
We claim that the following holds: 

(4.20) f AfD G jk {-,y)D a ^ = ^{y), Vc/> G (fl) N . 

To see (|4.20j) . write </> = ?70 + (1 — ??)</>, where rj G C 1 ^ (B dy (y)) is a cut-off function 
satisfying 77 = 1 on B dy/2 (y). Then, ((OJ), (|4~T7) . and P~T9"j) yield 

4> k {y)= lim + r/^ fc + lim / (1 - ry)0 fc 

= lim / A^D^^(-,y)Z> a (# i )+ lim / A^D P G% y)D a ((l 
= / ^/^ /3 G jfe (-,y)^ Q (^)+ / J^DpG jk {-,y)D a {{l-r,)^) 

= J AfD p G jk ( • , y)D a ^ as desired. 

Next, we claim that G( ■ , y) = on dfl in the sense that for all 77 G C£°(f2) satisfying 
77 = 1 on B r (y) for some r < d y , we have 

(l-^-^GYo 1 ' 2 ^)™ 
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To see this, it is enough to show that 

(4.21) (i- v )GP^.,v)^{l-Ti)G(-,y) in Y 1 - 2 (Sl) N * N as n -> oo, 

for (1 - rf)G p "{ ■ ,y) E Yq' 2 (Q) NxN for all fjt > 1 and Yq' 2 (Q) is weakly closed in 
Y 12 (n) by Mazur's theorem. To show ([4~2il . we note that (|4TT9]) yields 

(l-T7)G w (',y)^= / G kl (-,y)(l-r,)<f> = lim / ,y)(l - r,)cj> 

Jn Jn 



lim / (l-r))G%(-,y)ct>, e L*T*(fi), 



D((l-r))G kl (-,y))-il> = - / G u (- ,y)Dri^ + / DG W ( • , y) • (1 - t?)i/> 

■/o Jo 

= - lim / G%(-,y)D V -<4>+ Urn / DG£j* ( • , y) • (1 - r?) V 
v->oGj n v^°°Jn 

= lim / £>((l-»,)G#(-,tf)).tf, V^eL 2 (Of. 

By using the same duality argument as in (|3 . 29[) . we derive the following esti- 
mates that correspond to (|3.30|) - (|3.36p : 

(4.22) \\G(-,y)\\ LP{Br(y)) <C p r 2 ~ n+n /P, Vr < d yi Vp e [1, ^), 

(4.23) ^pr 1 - 1 ^, Vr<d y , Vpe[l,^), 
(4-24) l|G(-,y)|| rl , 2(nVBrW) <Cr x - n / 2 , Vr < d v /2, 

(4.25) |{ajen:|G(s,y)|>t}|<Ot - ^, Vt > (d y /2) 2 - n , 

(4.26) |{aGn:|D a .G(x,y)|>i}|<C7*-5*r, Vt > {d y /2)^ n . 

Also, we obtain pointwise bound and Holder continuity estimate for G( ■ , y) 
corresponding to (|3 . 38[) and (|3.39[) . respectively, as follows. Denote by v T the fc-th 
column of G( • , y) and set R := d x y /2, where 

(4.27) d Xt y := mm(d x ,d y , \x - y\). 

Since v is a weak solution of Lu = in B 3R / 2 (x) C O \ B R / 2 (y), it follows from 
([2~TU]) and P~2"Ij) that 

\v(x)\ < CR^' 2 \\v\\ L2 * ((ABR/2(y)) < CR 2 - n , 

which in turn implies that 

(4.28) \G(x,y)\ < Cd 2 ,'™, where d x>v := mm(e2 a , d y , \x - y\). 
In particular, we have 

(4.29) \G(x,y)\ <C\x~y\ 2 - n if \x - y\ < d x /2 or \x - y\ < d y /2. 
Similarly, it follows from (l2~TTj) and P~2"4]) that 

(4.30) [v] 2 cm{BR ( x )) < CR 2 - n - 2 ^> I \Dv\ 2 < CR 2 ^-^. 
Therefore, we find that 

(4.31) \G(x,y)-G(z,y)\<C\x-zrd 2 ~ y n -^ if \x - z\ < d x , y /2, 
where d XtV is given by (|4.27|) . 
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Denote by t G° ( • , x) the averaged Green's matrix of *L in f2 with a pole at x G CI. 
Observe that we have an identity corresponding to (|3.42|) . 

(4-32) / t G% l {-,x) = -l G p lk {-,y). 

Let *G( ■ , x) be a Green's matrix of *L in with a pole at x G f2 that is obtained 
by a sequence {ov}^^ tending to 0. Then, by a similar argument as appears in 
Section [3.6i we obtain 



(4.33) G lk {x,y) = t G kl {y,x), Vfc, I — 1, ... ,N, Vz.yefi, x ? y, 
which is equivalent to say 

(4.34) G(x,y)= t G(y,x) T , Vx,yGCl, x + y. 

Using (|4.34|) . we find that G(x, ■ ) satisfies the estimates corresponding to (|4. 22[) 
(|4~26f and (|4~3T|l . Moreover, by following the lines (|3~44| - ((3~45| and using (|4.32|) 
we obtain 

(4.35) G p (x,y) = ~f G{x,z)dz. 
Therefore, by the continuity, we find 

(4.36) lim G p (x,y) =G(x,y), VxjgH, x ? y. 

Finally, we summarize what we obtained so far in the following theorem. 

Theorem 4.1. Let CI be an open connected set in R ra . Denote d x := dist(a;, 9f2) 
for x G CI; we set d x = oo if CI = M™. Assume that operators L and L satisfy the 
property (H). Then, there exists a unique Green's matrix G(x,y) = (Gy (x, 
(x, y £ Ct,x 5^ y) which is continuous in {{x, y) G CI x Cl : x ^ y} and such that 
G(x, ■ ) is locally integrable in CI for all x G CI and that for all f = ( f 1 , . . . , f N ) T G 
C£°(Q) N , the function u = (it 1 , . . . , u N ) T given by 

(4.37) u{x) := / G(x,y)f(y)dy 
belongs to Yq' 2 (CI) n and satisfies Lu — f in the sense 

(4.38) / A&DpulDcP = f f l <f>\ V0 G C?(0) N . 
Jn Jn 

Moreover, G(x,y) has the properties that 

(4.39) f Afn G ]k {-,y)B a <\> 1 = <f> k (y), Vcb G C™(Cl) N 
Jn 

and that for all r\ G C^°(Cl) satisfying rj = 1 on B r (y) for some r < d y , 

(4.40) (l- V )G(-,y)EY 1 ' 2 (Cl) N * N . 
Furthermore, G(x, y) satisfies the following estimates: 

(4.41) \\G(-,y)\\ LP(Br(y)) <C p r 2 - n+ ^, Vr<d v , Vp G [1, 

(4.42) \\G(x, -)\\ LHBAx)) <C p r 2 - n+n l p , Vr < d x , Vp G [1, 

(4.43) pG ! (-,y)|| i p(B r c w ))<Cpr 1 - B+n / 1 ', Vr < d„ Vp G [1, ^), 

(4.44) ||I>G(x, Olli,^))^^^"^. Vr<4, Vp G [1, 
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(4-45) l|G(-,y)|| ri , Wr(y)) <Ct^I\ Vr<d y /2, 

(4-46) \\G(x, Olly^n^*)) < Or 1 -"/ 2 , Vr < d x /2, 

(4.47) \{xen: \G{x, y)\>t}\< Ct~^ , Vi > K/2) 2 "™, 

(4.48) |{y GQ: \G(x,y)\ > t}\ <Ct~^, Vt > (d x /2) 2 - n , 

(4.49) \{xen:\D x G(x,y)\>t}\<Cr^, Vt > (d y /2) 1 - n . 

(4.50) \{yefl:\D v G(x,y)\>t}\<Ct-^, Vt > (d x /2) 1 ~ n , 

(4.51) \G(x, y)\ < C(j?~y, where d XiV := mm(d x ,d y ,\x - y\), 

(4.52) \G(x,y)-G(z,y)\<C\x-zrd 2 - y n -^ if \x - z\ < d x , y /2, 

(4.53) \G(x,y)-G(x,z)\<C\y-zrd 2 x - n -»° if \y - z\ < d x , y /2, 

where C = C(n, N, A, A, /i , i? ) > and C p = C p (n, N, A, A, H ,p) > 0. 

Proof. Let G p (x, y) and G(cc, y) be constructed as above. We have already seen that 
G is continuous on {(x, y) 6 il x £1 : x ^ y} and satisfies all the properties (|4.39p 
- (|433l) . Also, as in the proof of TheoremO we find that for all / £ C™(fl) N , 
there is a unique u G (Fq 1 ' 2 ^) n C(fi)) Ar satisfying 

/ AffDpuWaV* = [ /V, V« G Fq 1 ' 2 ^)^. 
Jn Jn 

If we set = G p kl {x, • ) above, then by (43]) and (|4~34|) . we find 



Fix r < d x /2. By (|4"T7|) . gUJ), and (14TM)) . we have 



(4.54) / G2i(x, -)f = / ^A^-,*)^' = + « 



lim / G&(s, Of = lim / G#(*, -)f + / G&(x, Of 



G ki (x,-)f+ G ki (x,-)f 

Bt{x) if!\Bi(i) 

- / G M {x, Of- 



Therefore, (|4.37p follows by taking the limits in (|4.54|) . By proceeding as in the 
proof of Theorem 13. 11 we also derive the uniqueness of Green's matrix in f2. □ 

4.2. Boundary regularity. Let E be any subset of f2 and u be a W 1,2 {£1) function. 
Then we shall say u = on E (in the sense of W 1,2 (fl)) if u is a limit in VF 1,2 (f2) 
of a sequence of functions in C£°(f2 \ E). 

We shall denote ^r(x) := dil Pi Br(x) for any R > 0. We shall abbreviate 
= Qr(x) and E^ = T,r(x) if the point a; is well understood in the context. 

Lemma 4.2 (Boundary Poincare inequality). Assume that \Br > 9\Br\ for 
some 6 > 0. Then, for any u G W 1 ' 2 ^^) satisfying u = on E^, we have the 
following estimate: 

(4-55) IMIi=(n R ) < S^II^ILw 
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Proof. Since u = in we may extend u to a W 1,2 (Br) function by setting 
u = in S := 5r \ fi. Note that Du = in 5. Then the lemma follows from (7.45) 
in PH p. 164]. □ 



Lemma 4.3 (Boundary Caccioppoli inequality). Let the operator L satisfy condi- 
tions (|2.2I) . (|2.3|> . Suppose u is a W 1,2 (Qr) n solutions of Lu = in £Ir satisfying 
u = on Tir. Then, we have 

(4.56) \\Du\\ L , {nr) <^- r \\u\\ L2(QR) , V0<r<i?, 

where C = C(n,N,X,A) > 0. 

Proof. It is well known. □ 



Definition 4.4. We say that f2 satisfies the condition (S) at a point x £ dtt if 
there exist 9 > and i? a G (0, oo] such that 

(4.57) \B R (x)\n\>d\B R (x)\, VR<R a . 

We say that O satisfies the condition (S) uniformly on £ C <9S1 if there exist 9 > 
and i? a such that (I4.57|) holds for all x € S. 

Definition 4.5. Let fi satisfy the condition (S) at i 6 9f2. We shall say that 
an operator L satisfies the property (BH) if there exist Hi, Hi > such that if 
u £ W 1,2 (Sl]i(x)) N is a weak solution of the problem, Lu = in £Ir(x) and u = 
on E^(x), where R < R a , then u satisfies the following estimates: 

(4.58) / < ffi(-) / |L>«r, V0<r<s<R. 

Lemma 4.6. There exists eo = eo(n, Ao, Ao) > smc/i £/ia£ if the coefficients of the 
operator L in (|2.1[) satisfies (|2.8[) m Lemma \2.'A then L satisfies the property (BH) 
with [i\ — /ii(n, Ao, Ao, 6*) > and H\ = Hi(n, N, Ao, Ao, 9) > 0. 

Proof. Throughout the proof, we shall abbreviate f2 r = Q r (x) for any r > 0, 
S r = S r (a;), the point x £ dfl to be understood. For any s < R < R a , let 
v l (i — 1,...,N) be a unique W^ 1,2 (ri s ) solution of Lqv 1 — in il s satisfying 
v l - u l £ Wo' 2 (Sl s ), where L v i = -D a (a a/3 D /3 v i ). 

We claim that there exist ^{n, Ao, A , 9) > and C(n, Ao, A , 9) > such that 
the following estimate holds: 

(4.59) / \Dv\ 2 <C(-j / |L>v| 2 , V0<r<s. 

We first note that we may assume that r < s/8; otherwise (|4.59|) becomes trivial. 
Since each v l satisfies v l = on S s , it follows from Theorem 8.27 [TOl PP- 203-204] 
and Theorem 8.25 [TO; pp. 202-203] that there is fi 2 = ^{n, A , A , 9) > and 
C = C(n, Ao, A , 9) > such that 

(4.60) oscv 1 < Cr^s-^ sup |«*| < Cr' i2 s-^-"/ 2 ||^|| L 2 (f2 } . 
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In particular, the estimate (14.60|) implies v t {x) = lim v l (x) = 0. Then, Lemma l4~3l 

X — *x 

and Lemma [4721 imply that for alii = 1, . . . , N (recall r < s/8) 



\Dv*\' <Cr~ 2 I |«*| a = Cr~* \ \v l -v\x)\ 

2 



<>_, 



< Cr n - 2 { osc v*\ <C(- yl 2+2/18 ---^ / " ' ~ 



n 2T - / V s 



S / V 



and thus we have proved the claim. 

Next, note that w := u — v belongs to W ' (Q S ) N and thus it satisfies 

xf \Dw\ 2 < [ a af3 D fj w l D a w l = [ (a af3 5 l3 - Af)D p u 3 D a w\ 

Therefore, we have 

(4.61) / \Dw\ 2 <(X- x \\e\\ Loa f [ \Du\ 2 , 

where e(x) is as defined in (|2 . 8[) . By combining (|4.59|) and (|4.61[) . we obtain 
/ \Du\ 2 < C (-) n ~ 2+2 ^ [ \Du\ 2 + Cv\\ef L ^ [ \Du\ 2 , V0 < r < s. 

Now, choose a /ii S (O,/^)- Then, from a well known iteration argument (see, e.g., 
[H Lemma 2.1, p. 86]), it follows that there is eo such that if ||e|| ioo < eo, then 
(1438)1 holds. □ 

Theorem 4.7. Let the operator L satisfy the properties (H) and (BH). Assume 
that Q satisfies the condition (S) at x S dQ with parameters 8, R a . Let x G f2 such 
that \x — x\ = d x < R/2, where R < R a is given. Then, any weak solution u of 
Lu = in Qr(x) satisfying u — on ~E,n(x), we have 

(4.62) |u(a:)| < Cd^R 1 -^ 2 ^ \\Du\\ L2{nn{s)) , d x := dist(x, 9fi), 

where C = C[n, N, A, A, 9, /io, Hi, Hq, Hi) > and \i = min(/zo, /Ui). 

Proof. The proof is an adaptation of a technique due to Campanato 0]. In this 
proof, we shall use the notation u„ := £ n u. Also, we shall abbreviate d = d x . 
Observe that 

(4.63) n d (x) = B d {x) c n 2d (x) n n 2d {x). 

We may assume that R > 3d so that Sl 2d (x) C £Ir(x); otherwise 2d < R < 3d and 
(|4.62p follows from Lemma \2. 41 We estimate u(x) by 

\u(x)\ < \u(x) - u xM \ + \u xM - u xM \ + \u s ,2d\ := L + IL + LLL. 

We shall estimate / first. For any r\ < r 2 < 2d, we estimate 

(4.64) \u x , ri - u Xtr2 1 2 < 2 \u(z) - u x<ri \ 2 + 2 \u(z) - u x<T2 1 2 . 
Note that since B d (x) C f2, we have 

\Vl r {x)\>Cr n , Vr<2d. 
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Therefore, by integrating (|4.64[) over Q ri (x) with respect to z, we estimates 
(4.65) \u x , ri - u X: r2 \ 2 < Crf" / |it - u x , ri \ 2 + / \u - u x ^ | 2 • 



Since tt = on E#(a;), we may extend u to Bji(x) as a W ' function by setting 
m = 0oii Bfj (x) \ Q. Therefore, by a version of Poincare inequality (see, e.g. (7.45) 
in [TQl p. 164]), we have for all r < 2d, 

(4.66) I \u-u x , r \ 2 <[ \u ~ u x , r \ 2 < Cr 2 [ \Du\ 2 = Cr 2 [ \Du\ 2 . 
Ja r J B r JB r Jn r 

Therefore, by (|4.65|) and (|4. 66[> . we obtain 

2,2/ in„.|2 



(4.67) \u x . ri - u xr2 1 < Of" rf / |£>«| z + r z 2 / 

V isVj(j) Jn T2 (x) 

Next, we claim that the following estimate holds: 



(4.68) / \Du\ 2 < C [ - ) / \Du\ z , Vr < 2d 



r \ n — 2+2fi r „ 



RJ 



We first consider the case when r < d. Note that in this case, we have O r (x) = B r (x) 
and fld(x) — Bd(x). Since L satisfies (H), it follows from (|4.63|) that 

(4.69) / \Du\ 2 <C(-) / \Du\ 2 <C[-) / \Du\ 2 . 
On the other hand, since L satisfies (BH), it follows from (|4.58[) that 

(4.70) / |^| 2 <C - / |i?«| 2 . 

By combining (|4.69|) and (|4.70|) . we obtain (|4.68[) . Next, consider the case when 
d < r. In this case, we have £l r (x) C ^2r{x), and thus it follows from (|4.58|) 



9 f ,9 / r \n-2+2;i 

Ditr < / LDtt < CI - / LDu 



2 



We proved the claim (14. 68[) . 

Now, by using (|4. 68[) . we estimates (|4.67p as follows (recall r\ < r 2 < 2d): 

(4.71) K, ri - u x , r2 \ 2 < Cr^ n (r^ + r^R 2 -^ f \Du\ 2 . 

Jn R (x) 

For any r < 2d, set r\ = r2~( !+1 ) and T2 = r2~ l in (14.7ip to get 

\u x r2 _ (i+1) - U Xir2 -i\ 2 < Cr 2 ^2- 2 ^R 2 - n - 2 ^ [ \Du\ 2 . 

Jn R (x) 

Therefore, for < j < k, we obtain 

k-l 

\ u x,r2- k ~ u x,r2-i \ — ^ ' | u x,r2~ — u x,r2~ i \ 

i=3 

(4 72) l °° \ 



C2 -3K r K R l-n/2-n ,, D < 



L 2 (n R (x)) 
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By setting r — 2d, j — 0, and letting k — > oo in (|4.72|) . we obtain 

(4.73) I = \u(x) - u xM \ < Cd^R 1 -^ 2 ^ \\Du\\ L2{nR{s)) . 

Next, we estimate III. Since \B r (x) n Bd{x)\ > Cr n for r < 2d, we have 

(4.74) |fir(*)| > Cr", Vr<2d. 

Also, as in (|4. 66[> . we have for all r < 2d (recall u = on Br(x) \ ft) 

(4.75) / \u-u £ . r \ 2 <{ \u - u s , r \ 2 < Cr 2 [ \Du\ 2 = Cr 2 [ \Du\ 2 . 



Jn r J B r J B r Jn r 

Therefore, as in (|4.67|) we have for n < r 2 < 2<i, 

|us, n -« S ,, 2 | 2 < Cr- n Irf f \Du\ 2 + r 2 2 f \Du\ 2 ) . 

y Jn ri (x) Jn r2 (x) J 

Then, by using the property (BH), we obtain (c.f. ljl7T2"]) . (|4TF5]) ) 

(4.76) - u sM \ < Cd^R 1 -^ \\Du\\ L2(nRm , 

where u(x) := lim^oc u s 2-fc r - (note that (|4.72[) implies exists). It follows 

from P~fg]) . (|4"35)) . and (|4"35)) that for any r < 2d, we have 

Iw,.,| 2 < / \u\ 2 < Cr- n [ \u\ 2 



Q r (x) JQ r (x) 



< Cr 2 ~ n 



[ \Du\ 2 <Cr 2 - n (^ 2+2 " [ \Du\ 
Jn T {x) KK/ Jn R (s) 



= Cr 2 ^R 2 - n - 2 » / \Du\ 2 , 

and thus that u(x) — 0. Therefore, by (|4.76|) we obtain 

(4.77) III = \u sM \ = \u(x) - u mfi d\ < Cd^R 1 -^ 2 -^ \\Du\\ L2inR{s)) . 
Finally, we estimate II. 

(4.78) \u X: 2d ~ Ux,2d\ 2 < 2 \ U ( Z ) - u xad\ 2 + 2 \ U ( Z ) ~ u x,2d\ 2 ■ 

By integrating (|4.78|) over B d (x) C f^dM n ^ 2 d(aO with respect to z, we estimate 

k x ,2d - Mx,2d| 2 < CgT" / \U - U xa d\ 2 + / |«-«x,2d| 2 



(4-79) <Cd 2 - n (f \Du\ 2 +f \Du\ 2 ) 

< Cd 2 ^R 2 - n - 2 ^ [ \Du\ 2 , 

where we have used ()4.66|) . (|4.75|) . (|4.68[) . and (I4.58| . Therefore, by combining 
(|473l) . (I477j) . and (I4J9]) . we obtain (|4T62j) . □ 

Theorem 4.8. Let t/ie operators L, L satisfy the properties (H) and (BH). Assume 
that satisfies the condition (S) uniformly on <9f2 with parameters 0,R a . Denote 

R x , y := min(|a: - y \ , 4i? a ). 
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Then the Green matrix G(x,y) satisfies 

(4.80) \G(x,y)\ < Cd^Ri^l^d 1 -^ 2 if d x < R x , y /8, 

(4.81) \G(x,y)\<Cdl£R£?' 2 -'>dl-»'* ifd y <R x , v /8, 

where C — C(n, N, A, A, 9, /i , Hi, H , Hi) > and /j = min(/j ,Mi)- As a conse- 
quence, we have G( ■ , y) = 0, G(x, ■ ) = on dVL in the usual sense. 

Proof. We only need to prove P~8"0"]) . for P~gTj) will then follow from (1434]) . Set 

i? = R x ,y/4, r = dy/2, and choose a; £ <9f2 such that |x — x\ = d x . Then, since 

d y < \x - y\ + d x < | \x-y\ , 

we have 

\y - x\ > \x - y\ - d x > | \x - y\ > R + r, 

and thus, iln(x) C ft\B r (y). Now, we apply Theorem l4.7l with u = G(- ,y). Then, 
by (|432J) and fl~2"4")) , we obtain 

\G(x,y)\ < Gd^R 1 -^ \\DG(-,y)\\ L2{nXBr{y)) < Cd^E^^d^^ . 

The proof is complete. □ 

Remark 4.9. We note that in the scalar case, the maximum principle yields (see 
[Til Theorem 1.1]) 

(4.82) G(x,y) < C \x - y\ 2 ~ n , Vx ^ y £ fi. 

Then, by the boundary Caccioppoli inequality, we have (c.f. (|4.4|) - (|4.7|l ) 

\DG(-,y)\ 2 < Cr 2 ~ n , Vr > 0. 



/ 



n\B r (y) 

Therefore, in the scalar case we don't need to require that r < d y /2 (or r < d x /2) 
in the proof of Theorem 14.81 and we may as well set r = \x — y \ /2 to get 

G(x, y) < CdZRl; y n / 2 ~» \x - y^ 2 if d x < R x , y /8, 

G(x, y) < Cd^R 1 -™/ 2 -* \x - y\ l - n/2 if d y < R x>y /8. 

In particular, if G(x, y) is the Green's function on M™ , then we obtain 

G(x, y) < Cd% \x - y\ 2 - n -» if d x < \x - y\ /8, 
G(x, y) < Cd» \x - y\ 2 - n -<* if d y < \x - y\ /8, 
for <9R™ satisfies the condition (S) with 9 = 1/2 and R a = oo. 

5. Remarks on VMO coefficients case 

Definition 5.1 (Sarason [IH])- For a measurable function / defined on R™, we 
shall denote f xr = f B , , f and for < S < oo we define 

(5.1) M s (f) ■■= sup sup / |/ - 7 Xir | ; M (f) := lim M s (f). 

r<S JB r ( X ) 

We shall say that / belongs to VMO if M (f) = 0. 
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Definition 5.2. We say that the operator L satisfies the property (H) ;oc if there 
exist no, Ho, R c > such that all weak solutions u of Lu = in Br = Br(xq) with 
R < R c satisfy 

f o / r \n-2+2fj, r 

(5.2) J \Du\ 2 <H Q (-) J \Duf , 0<r<s<R. 

Similarly, we say that the transpose operator *L satisfies the property (H) /oc if 
corresponding estimates hold for all weak solutions u of l Lu = in Br with R < R c . 

Lemma 5.3. Let the coefficients of the operator L in (|2.1[) satisfy the conditions 
(|2.2p and (|2.3[) . If the coefficients belong to VMO in addition, then L satisfies the 
property (H) Zoc . 

Proof. It is well known that if the coefficients are uniformly continuous, then L 
satisfies the property (H) Zoc ; see e.g. [HI pp. 87-89]. Essentially, the same proof 
carries over to the VMO coefficients case. One only needs to make a note of the 
following two facts. First, a theorem of Meyers [16] implies that there is some 
p = p(n, N, A, A) > 2 such that if u is a weak solution of Lu = in Br(x), then 

l/p , v 1/2 

< C I I \Duf , Vr < R/2. 
Secondly, note that the John-Nirenberg theorem [T3] implies that 

(7 |/-7r,*rl <C(n,g)M s (f), Vr<c(n)d, V 9 e(0,oo), 

\JB r (x) J 

where M${f) is defined as in (|5.ip . For the details, we refer to [3j pp. 47-48]. □ 

In the rest of this section, we shall assume that the operators L and *L satisfy 
the property (H) ioc with parameters /xq, Ho, R c - We shall denote 

(5.3) r x := min(d x , R c ), f x ^:=miri(d x , y ,R c ), 

where d x = dist(x, dfl) and d XtV is as in (|4.28|) . It is routine to check that all 
estimates appearing in Section [4.11 remain valid if d x , d x%v are replaced by r x , f x ^ y , 
respectively. Therefore, we have the following theorem: 

Theorem 5.4. Let Q be an open connected set in W l . Denote d x := dist(a;, dQ) for 
x G rt; we set d x = oo if fl = W 1 . Assume that operators L and *L satisfy the prop- 
erty (H) (oc . Then, there exists a unique Green's matrix G(x,y) — (Gi -(x, y))f^- =1 
(x,y G Ct,x 7^ y) which is continuous in {(x,y) G £1 x ft : x ^ y} and such that 
G(x, ■ ) is locally integrable in f2 for all x G f2 and that for all f = (f 1 , . . . , f N ) T G 
C^°(tt) N , the function u = (it 1 , . . . , u N ) T given by 




(5.4) u{x) := / G{x,y)f(y)dy 
belongs to Y 1 ' 2 (ri) A ' and satisfies Lu — f in the sense 

(5.5) / AffDpuW^ = f f<f>\ V0 G C?(n) N . 
Jn Jn 

Moreover, G(x,y) has the properties that 

(5.6) J^AffD f3 G 3k (-,y)D a ^=^ k (y), V<t>£C™(Q) N 
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and that for all r\ 6 C^ c '(Sl) satisfying rj = 1 on B r (y) for some r < d y , 
(5.7) {l-r,)G(',y)€Y^(n) N ^. 

Furthermore, G(x,y) satisfies the following estimates: For r Xl r y ,f x ^ y as in (|5.3j) . 

5.8) \\G(-,y)\\ LP{BAy)} <C p r 2 - n+n lv, Vr<r„ Vp G [1, 

5.9) IIGC^OIIlW,))^^^-" 4 *^, Vr<r X) Vp G [1, ^j), 

5.10) ||£>G(-,y)|| iP(flr(j/)) fCCpr 1 -'^, Vr < r,, Vp G [1, 

5.11) ||i?G( a : ) -)ll L -(fl rW )<Cpr 1 -' l+ »^, Vr < r Xl Vp G [1, 

5-12) l|G(-,tf)|| y i. a(n \ Br(s)) <Cr 1 -' l / 2 ) Vr < r„/2, 

5-13) |[G(a:, • )|| yl , Wr( ,)) < Cr 1 -"/ 2 , Vr < r x /2, 

5.14) \{x€Sl:\G(x,y)\>t}\<Cr^, Vt > (r y /2) 2 - n , 

5.15) \{y& SI :\G(x,y)\>t}\ <Ct-£*, Vt > (r x /2) 2 - n , 

5.16) \D x G(x, y)\ > t}\ < Ct~^ , Vt > (r^/2) 1 -™, 

5.17) \{yeSl:\D y G(x,y)\>t}\<Ct-^, Vt > (r^) 1 "™, 

5.18) <Cf 2 -", Vx,yGQ, 

5.19) |G(z,y)-G(*,2/)| <C\x-z\*> f 2 x ~ n -^ if \x - z\ < r B ,„/2, 

5.20) |G(a:,»)-G(a:,z)| <C|y-zp f 2 ""'^ i/ | y - z| < f^/2, 
w/iere C = C(n, N, A, A, fi , Hq) > and C p = C p (n, N, A, A, /j, , H ,p) > 0. 
Remark 5.5. Dolzmann-Miiller [B] derived a global estimate 

(5.21) < C|x-y| 2 -" Vx,yeSl, x + y, 

assuming that SI is a bounded C 1 domain. We have not attempted to derive the 
corresponding estimate here. However, we would like to point out that the constant 
C in their estimate depends on the domain (e.g., the diameter of the domain and 
also some characteristics of dSl) while our interior estimate (|5 . 18[) does not. 
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